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In re Application' of: 




L AMENDMENTS 



A. In the Specification 

Please enter rewritten paragraph 1 1 on page 3 to read as follows: 

Figure 3 shows a schematic drawing of a regulatory network associated with a reaction 
in a metabolic network. Integration of a stoichiometric model and a logical model is achieved 
through regulatory restraints (logic values of reaction processes) which are used to refine 



appropriate reaction constraints in the model. If rxnLogic = 1 then use Activity constraints; If 
rxnLogic 0 then use Inactivity constramts. Activity constraints set for rxnsioich: (lower bound 
= 0, upper bound = INF or #); Inactivity constraints for rxn stoich^ (lower bound = 0, upper 
bound = 0). Logic functions: ai = (activator/inhibtor) • TF; az = 1; C| = TF* • pri • genel • 

gene2; C2 = prj • geneS; 1 1 = Mgcnel; h = Mgene2; h = Mgcne3; pl = Pgenel ; • Pgene2 ' Pgene3; rSHugic 

= Protein • Cofactor • Substratei • Substrate:. Time delays can be specified for the switching of 
each memorization variable after a triggering change in the associated function. 



Please enter rewritten paragraph 15 on page 4 to read as follows: 

Figure 7 shows a schematic drawing of a simplified core metabolic network. Table 4 
provides the stoichiometry of the 20 metabolic reactions included in the network. 



Figure 8 shows, in Panel A, a simulation of aerobic growth of E. coli on acetate with 
glucose reutilization; and in Panel B, in silico arrays showing the up- or down-regulation of 
selected genes, or activity of regulatory proteins, in the regulatory network. Panel A shows 
three time plots showing experimental data (closed squares or triangles) and the corresponding 
simulations performed using the combined regulatory/metabolic model (solid lines) as well as 
the stand-alone metabolic model (dashed lines). Table 5 provides parameters used to generate 
the plots in Panel A. 
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•lease enter rewritten paragraph 16 on page 4 to read as follows; 




9^ 



±OOS.7 W'l, „iy'rj.r^'^4^- 
In re Applicatioaof: PATENT 
PalssonetaL AttoraeJIPcket No,: UCSD1330-2 

Application No.: 10/087,441 
Filed: March 1,2002 
Page 3 



please enter rewritten paragraph 17 on page 4 to read as foUowsT^ 
Figure 9 shows, in Panel A, a simulation of anaerobic growth of £. coli on glucose; and 
in Panel B, in silico arrays showing the up- or down-regulation of selected genes, or activity of 
regulatory proteins, in the regulatory network. Table 6 provides the parameters used to 
generate the plots in Panel A. 

^Please enter rewritten paragraph 18 on page 4 to read as followsP) 
Figure 10 shows, in Panel A, a simulation of aerobic growth of E. coli on glucose and 
lactose; and in Panel B, in silico arrays showing the up- or down-regulation of selected genes, 
or activity of regulatory proteins, in the regulatory network. Table 7 provides the parameters 
used to generate the plots of Panel A, Panel A shows time plots showing experimental data 
(triangles) and the corresponding simulations performed using the combined 
regulatory/metabolic model (thick solid lines), the stand-alone metabolic model (dashed lines), 
and the kinetic model described in Kremling (Metabolic Eng. 3:362-379 (2001)) (thin solid 
line). 



Please enter rewritten paragraph 139 on page 48 to read as follows: 



A skeleton of the biochemical reaction network of core metabolism was formulated, 
including 20 reactions, 7 of which are regulated as shown in Figure 7. This network provided a 
simplified representation of core metabolic processes including glycolysis, the pentose 
phosphate pathway, TCA cycle, fermentation pathways, amino acid biosynthesis and cell 
growth, along with corresponding regulation pathways including catabolite repression, 
aerobic/anaerobic regulation, amino acid biosynthesis regulation and carbon storage regulation. 
The skeleton biochemical reaction network was represented as a skeleton combined 
regulatory/metabolic model where reactions were represented as linear equations of reactants 
and stoichiometric coefiBcients and regulation was represented by regulatory logic statements 
as shown in Table 4. As shown in Figure 7 and Table 4, four regulatory proteins (Rpo2, 
RPcl,RPh and RPb) regulated 7 of the 20 reactions in the skeletal network and model. 
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Please enter rewritten paragraph 156 on page S3 to read as follows: 



E. coli has been observed in vivo to secrete acetate when grown aerobically on glucose 
in batch cultures; when glucose is depleted from the environment, the acetate is then reutilized 
as a substrate. Using the combined regulatory/metabolic and stand-alone metabolic models, 
activity of an aerobic batch culture of E. coli on glucose minimal medium was simulated. 
Panel A of Figure 8 shows three time plots showing experimental data (closed squares) and the 
corresponding simulations performed using the combined regulatory/metabolic model (solid 
lines) as well as the stand-alone metabolic model (dashed lines). In the acetate plot, the 
regulatory/metabolic model predictions differed from that of the stand-alone metabolic model, 
as shown. Table 5 provides the parameters required to generate the time plots where 
parameters were estimated or obtained from Varma and Palsson Appl. Env. Micro. 60:3724- 
3731 (1994). The major difference between the combined regulatory/metabolic and metabolic 
stand-alone simulations is in the delayed reaction of the system to depletion of glucose in the 
growth medium. The stand-alone metabolic networic is unable to account for the delays 
associated with protein synthesis. 



f5 



Please enter rewritten paragraph 159 on page 54 to read as follows: 



The in silico models were used to simulate anaerobic growth on glucose, the results of 
which are shown in Figure 9 which was generated using the parameters provided in Table 6. 
Under these conditions, the stand-alone metabolic model made similar predictions as the 
combined regulatory/metabolic model, with a notable exception: the combined 
regulatory/metabolic model was able to make predictions about the use of a particular isozyme. 
For example, both models require fumarase activity as part of the optimal flux distribution; 
however, of the two models only the combined regulatory/metabolic model was able to 
specifically determine that the fumB gene product [which as being] is expressed under 
anaerobic conditions. 
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^^ase enter rewritten paragr^h 160 on page 55 to read as follo^ ^T^ 
Aerobic growth of E, coli on glucose and lactose was simulated using the in silico 
models and compared to in vivo observations from mixed batch cultures and to results reported 
for a kinetic model as described in Kremling et al., Metabolic Eng. 3:362-379 (2001). Overall, 
the combined regulatory/metabolic model predictions were in good agreement with the in vivo 
observations, comparable with the predictions made by the Kremling model, and better than 
the predictions of the stand-alone metabolic model as shown in Figure 10 which was generated 
using the parameters provided in Table 7. The deficiencies in the ability of the stand-alone 
metabolic model to accurately predict the results of this experiment is most likely due to the 
concurrent uptake of glucose and lactose, resulting in much more rapid depletion of the 
substrates and a higher growth rate. Interestingly, because of the larger flux of carbon source 
uptake, the stand-alone metabolic model predicted that E. coli growth should be oxygen-, 
rather than carbon-limited in this case. Accordingly, the secretion of acetate and formate was 
predicted by the stand-alone metabolic model. In contrast, the combined regulatory/metabolic 
model predicted that no secretion will occur under these conditions. 

jt^ase enter rewritten paragraph 161 on page 55 to read as foUowsT]]) 
The in silico arrays for the sunulation (Figure lOB) showed one shift in gene 
expression, occurring just under five hours. The up-regulation of the lactose uptake and 
degradation machinery, together with key enzymes in galactose metabolism, enables the 
system to use lactose as a carbon source once the glucose in the medium has been depleted. 



Grayaiy\GP6294792.1 
101668-66 



In re Application of: 
Palsson et al. 



Application No.: 10/087,441 
Filed: March 1,2002 
Page 6 



J|| PATENT 

AttomejBRcketNo,: UCSD1330-2 



Please add the following: 

After the second page of Table 3 on page 62, please add the following tables: 



Table 4 provides the stoichiometry of the 20 metabolic reactions include in the network of 

Figure 7. 



REACTION 


name: 


REGULATION 


Metabolic Reactions 






-lA'lATP+lB 


Rl 




-1B+2ATP+2NADH+1C 


R2a 


IF NOT (RPb) 


-1C-2ATP-2NADH+1B 


R2b 




-IB+IF 


R3 




-IC+IG 


R4 




-1G+0.8C+2NADH 


R5a 


IFNOT(RPo2) 


-1 G +0.8 C +2 NADH 


R5b 


IFRPo2 


-I C +2 ATP +3 D 


R6 




-IC^NADH +3E 


R7 


IF NOT (RPb) 


-1G-1ATP-2NADH+1H 


R8a 


IF NOT (RPh) 


+1G+1ATP+2NADH-1H 


R8b 




-1NADH-1.02+1ATP 


Rres 


IFNOT(RPo2) 


Transport Processes 






-1 Carbonl +1 A 


Tel 




-1 Carbon2+l A 


Tc2 


IF NOT (RPcl) 


-1 Fext+l F 


Tf 




-1 D +1 Dext 


Td 




-lE+1 Eext 


Tc 




-IHext+lH 


Th 




-1 Oxygen +1 02 


To2 




Maintenance/Growth Processes 






-1 C-1 F-1 H-lOATP+1 Biomass 


Growth 




Regulatory Proteins 








RPo2 


IF NOT (Oxygen) 




RPcl 


IF Carbonl 




RPh 


IFTh 




RPb 


IFR2b 



GrayCaiy\G'T>6294792.l 
101668-66 



In re Application, of: PATENT 
Palssonetal. Attome^RketNo/. UCSD 1330-2 

Application No.: 10/087,441 
Filed: March 1,2002 
Page? 



Table 5 provides parameters used to generate the plots in Figure 8A. 



Parameters 






Initial Conditions 






[Biomass]o (g/L) 


0.003 


Estimated 


[Glucosejo (mM) 


10.4 


Estimated 


[Acetate]o (mM) 


03 


Estimated 


Strain-specific parameters 






Protein synthesis/degradation delay 
(hrs) 


0.4 


Estimated 


Biomass scaling factor 


1,3 


Varma and Palsson, 1994 


Uptake rate constraints (mlA/(gDCW*hr)) 






Glucose 


10.5 


Varma and Palsson, 1994 


Acetate 


2.5 


Estimated 


Oxygen 


15.0 


Varma and Palsson, 1994 


Table 6 provides the parameters used to generate figure 9A. 


Parameters 






Initial Conditions 






[Biomass ]o (g/L) 


0.002 


Estimated 


[Glucose]o (mM) 


10.5 


Estimated 


Strain-specific parameters 






Protein syntfaesis/degradation delay 
(hrs) 


0.4 


Estimated 


Biomass scaling factor 


1.3 


Varma and Palsson, 1994 


Uptake rate constraints (mM/(gOCW*kr)) 






Glucose 


18.5 


Varma and Palsson, 1994 


Oxygen 


0.0 


Varma and Palsson, 1994 
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Table 7 provides the parameters used to generate the plots of Panel lOA. 



Parameters 






Initial Conditions 






[Biomassjo (g/L) 


O.OU 


Estimated 


[Glucoselo (mM) 


1.6 


Estimated 


[Lactose]o (mM) 


5.8 


Estimated 


Strain-specific parameters 






Protein synthesis/degradation delay 
(hrs) 


0.5 


Estimated 


Unconstraining llacZ] (mmol/gDW) 


0.0274 


Estimated 


Constraining [lacZ] (mmol/gDW) 


0.0015 


Estimated 


Uptake rate constraints (mM/(gDCW*hr)) 






Glucose 


6.5 


Estimated 


Acetate 


3.0 


Estimated 


Oxygen 


15.0 


Estimated 
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